QED calculation of the 2pi/2-2s and 2p3/2-2s transition energies 
and the ground-state hyperfine splitting in lithiumlike scandium 

Y. S. Kozhedub\ D. A. Glazov\ A. N. Artemyev\ N. S. Oreshkina\ 
V. M. Shabaev\ I. I. Tupitsyn\ A. V. Volotka^ and G. Plunien^ 
^ Department of Physics, 
St. Petersburg State University, 
Oulianovskaya 1, Petrodvorets, 
St. Petersburg 198504, Russia 
^ Institut fUr Theoretische Physik, 
TU Dresden, Mommsenstrafie 13, 
D-01062 Dresden, Germany 

Abstract 

We present the most accurate up-to-date theoretical values of the 2pi/2''^^ ^^id 2p3/2-2s transition en- 
ergies and the ground-state hyperfine splitting in Sc^^+. All two- and three-electron contributions to the 
energy values up to the two-photon level are treated in the framework of bound-state QED without aZ- 
expansion. The interelectronic interaction beyond the two-photon level is taken into account by means of 
the large-scale configuration-interaction Dirac-Fock-Sturm (CI-DFS) method. The relativistic recoil correc- 
tion is calculated with many-electron wave functions in order to take into account the electron-correlation 
effect. The accuracy of the transition energy values is improved by a factor of five compared to the previous 
calculations. The CI-DFS calculation of interelectronic-interaction effects and the evaluation of the QED 
correction in an effective screening potential provide significant improvement for the 2s hyperfine splitting. 
The results obtained are in good agreement with recently published experimental data. 

PACS numbers: 12.20.Ds, 31.30.Jv, 31.10.+Z, 31.30.Gs 
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I. INTRODUCTION 



The dielectronic recombination process has proven to be a useful tool in high-precision mea- 
surements of the excitation energy of low-lying levels in middle-Z lithiumlike systems llli|2|]. By 
this method the energy of the 2p3/2-2s transition in Sc^^"^ was determined to be 44.3107(19) eV [2]. 
A significant improvement of the accuracy was announced recently by M. Lestinsky et al. ISi |4|], 
with the preliminary value of 44.3096(4) eV, and the work on further improvement of this value 
is in progress [iSJ. In these experiments the energy of the Rydberg resonances E^^s was measured. 
The Rydberg state energy -Ebind was evaluated by means of relativistic many-body perturbation 
theory (RMBPT). Then the excitation energy of the ion was determined as E^^c = -E'res + -E'bind- 
InRef. lli the theoretical value of i^exc for both 2pi/2-2s and 2p3/2-2s transitions was obtained 
by means of RMBPT, while for the quantum electrodynamic (QED) correction the result of Ref. 
^ was taken into account. The energy resolution achieved in these experiments also allowed for 
resolving the 2s hyperfine structure. As a result, the 2s hyperfine splitting of lithiumlike scandium 
was measured to be 6.21(20) meV 

The main goal of the present investigation is to evaluate the 2pi/2-2,s and 2p3/2-2s transition 
energies and the ground-state hyperfine splitting in lithiumlike scandium to the utmost accuracy 
aiming at a stringent test of the present state-of-the-art theoretical description of many-electron 
effects. Various contributions to the energy of the 2p-2s transitions are considered in the next Sec- 
tion. In order to meet the experimental accuracy, rigorous quantum electrodynamic calculations of 
the first two orders of perturbation theory are combined with large-scale configuration-interaction 
Dirac-Fock-Sturm (CI-DFS) calculations of the third- and higher-order contributions within the 
Breit approximation. The relativistic nuclear recoil corrections are calculated as well. The evalua- 
tion of the hyperfine splitting is accomplished in Section |llll The CI-DFS method is employed to 
obtain correlation effects of order l/Z"^ and higher. The radiative correction to hyperfine splitting 
is calculated with an effective local screening potential. 

Relativistic units are used throughout the paper {h = c = 1). 

n. 2j9i/2-2s AND 2j93/2-2s TRANSITION ENERGIES 

We start with the Furry picture, where in the zeroth-order approximation noninteracting elec- 
trons are bound by the Coulomb field of the nucleus. The Dirac equation yields zeroth-order 



2 



energies of the one-electron states. The homogeneously-charged-sphere model of the nucleus is 
employed with the value of rms radius (r^)^/^ = 3.5443(23) fm [|7|]. 

In leading order of the perturbation theory, diagrams of self-energy, vacuum polarization, and 
one-photon exchange arise. Techniques for the evaluation of these corrections nonperturbative 
in oiZ have been described in numerous publications (see, e.g., Ref. [SD). For the self-energy 
correction we inter pol ate the values presented in Ref. |^ for the 2 s and 2pi/2 states and those 
presented in Ref. lIlOll for the 2^3/2 state. The vacuum-polarization and one-photon exchange 
corrections are recalculated in the present work with inclusion of finite-nuclear-size effects. 

The second-order contributions can be classified as one-electron two-loop QED corrections, 
two-electron QED corrections, and two-photon exchange. Rigorous calculation of all two-loop 
QED corrections is a challenging problem. To date, the dominant part of these corrections was 
calculated in a wide range oi Z = 10 — 92 for the Is state only (see Ref. [11] and references 
therein). Recently, the corresponding results for 2s, 2]9i/2 and 2^3/2 states were presented for 
high-Z ions [12]. However, since for low values of Z the numerical evaluation of the second- 
order self-energy correction becomes rather difficult, so far one has to rely on the aZ expansion, 
which reads 
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where L = ln[(aZ)~^]. The values of the coefficients for the 2s state can be found in Appendix A 
of Ref. 1I13I1 and for the 2pi/2 and 2^3/2 states in Ref. [14]. Since the convergence of the expansion 
in aZ is known to be rather bad, we assume the uncertainty to be about 50% in our case. 

The two-electron QED corrections are represented by the diagrams of the screened self-energy 
and the screened vacuum-polarization. Rigorous evaluation of the screened self-energy in Li-like 
ions was performed in Ref. [15] for 2s and 2pi/2 states and in Ref. [Ilql for the 2^3/2 state. The 
screened vacuum-polarization correction was calculated in Ref. lilTll . We obtain the corresponding 
values for Z = 21 employing the procedure presented in these works. In order to estimate higher- 
order (in 1 /Z) terms of the screened QED correction, the following approximate scheme is used. 
The first-order QED correction is evaluated in an effective screening potential and the higher- 
order terms are extracted by subtracting the zeroth- and first-order terms. The uncertainty of the 
higher-order screened QED correciton obtained in this way is assumed to be 100%. 

The two-photon exchange correction is evaluated within the framework of QED, following our 



previous investigations IIS 

In order to evaluate the interelectronic-interaction corrections of third and higher orders we 
proceed as follows. The Dirac-Coulomb-Breit equation within the no-pair approximation is solved 
by means of the large-scale CI-DFS method 11201 12111 yielding the many-electron wave functions 
and the energy values. The interelectronic-interaction operator employed in the Dirac-Coulomb- 
Breit equation reads 



Breit 
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where a scaling parameter A is introduced in order to separate terms of different order in 1/Z 
from the numerical results with different A. Here i,j enumerate the electrons and o; is a vector 
incorporating the Dirac matrices. In this way, for small A, the total energy of the system can be 
expanded in powers of A, 



E{X) =Eo + EiX + E^X^ + J2 ^fcA' 
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The higher-order contribution E^^ = Xlfcla is calculated as E^^ = E{X = 1) — Eq — Ei — E2, 
where the low-order terms Eq, Ei, and E2 are determined numerically according to Eq. dH). Com- 
parison of El and E2 with corresponding QED results allows us to conclude that the uncertainty 
of the higher-order contribution due to the Breit approximation is less than 0.1%. 

The full relativistic theory of the nuclear recoil effect can be formulated only in the framework 
of QED [22]. To evaluate the recoil effect within the lowest-order relativistic approximation one 
can use the operator (see, e.g., Ref. Q) 
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where M is the nuclear mass and Pi is the momentum operator acting on the ith electron. The 
expectation value of Hm on the many-electron wave function of the system, obtained by the CI- 
DFS method, yields the recoil correction to the energy levels in all orders oi \/Z within the 
{aZYm? / M approximation. The electron-correlation effects contribute to about 20% of the total 
value and have to be taken into account in order to achieve the desirable accuracy. The one- and 
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two-electron recoil corrections of higher orders in aZ are taken from Refs. 11231 12411 . The recoil 
correction of the next order in m/M is negligible in the case under consideration. 

All contributions to the transition energies considered above are collected in Table I. For com- 
parison, previous theoretical results and available experimental data are presented as well. As one 
can see from the table, the theoretical values of the transition energies reported in this paper are 
about five times more precise than those in Ref. [120 and agree well with the experiments. Fur- 
ther improvement of the theoretical accuracy can be achieved by more accurate calculations of the 
higher-order screened QED effects. 



ni. HYPERFINE SPLITTING 



The ground- state hyperfine splitting of a lithiumlike ion is conveniently written as 

1 
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where nip is the mass of the proton, ji and / are the nuclear magnetic moment and spin, and fi^ de- 
notes the nuclear magneton. The one-electron relativistic factor A{aZ) can easily be derived from 
the Dirac equation utilizing virial relations [|26|l . The finite-nuclear-size correction 5 is evaluated 
numerically employing the homogeneously-charged-sphere model for the nuclear-charge distri- 
bution. The Bohr-Weisskopf correction e, arising due to the nonpointlike nuclear magnetization 



distribution, is evaluated within the single-particle nuclear model [|27 



m. 



The first-order interelectronic-interaction correction described by the function B(aZ) is eval- 
uated in the rigorous QED approach [|29|l . The dual-kinetic -balance (DKB) approach [|30(| is em- 
ployed to construct the complete set of one-electron wave functions from the B splines. The finite 
distributions of the nuclear charge and the nuclear magnetization are taken into account. The latter 
is introduced via the replacement of with F(r)/r^ in the hyperfine interaction matrix ele- 
ments. The explicit form of the function F{r) can be found in Refs. Oil 13211 . The higher-order 
correction C{Z, aZ)/Z'^ is obtained in the framework of the large-scale CI-DFS method. 

The QED correction Xrad is evaluated in one-loop approximation with an effecitve non- 
Coulomb binding potential Ves, which partly takes into account the interelectronic-interaction 
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effects. It is taken in the following form [|33 

Ves(r) = Kuc(r) + a 
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Here p is the total electron density, including the (Is)^ shell and the 2s electron. The parameter 
Xq, is taken to be a;^ = 2/3, which corresponds to the Kohn-Sham potential. To provide a proper 
asymptotic behavior, the potential V^s should be corrected at large r [35]. The one-electron spec- 
trum of the Dirac equation with V^g is constructed by means of the DKB method J3(j]. Since the 
potential V^g is assumed to be self-consistent, the standard iteration procedure is employed. The 
calculations performed are very similar to our recent calculations of the one-loop QED corrections 
to the g factor of Li-like ions [|36(1 . We mention also that the evaluation of the QED corrections to 
the hyperfine structure with an effective screening potential was performed in the past for the case 
of lithiumlike bismuth OTN . 

The individual contributions to the hyperfine splitting in litiumlike scandium are listed in Table 
II. For each contribution the corresponding term in the square brackets in Eq. ^ is explicitly 
written. For comparison, the experimental value from Ref. ||4|] as well as the previously published 



results by Shabaev et al. [13811 and by Boucard and Indelicato [|39h are presented. The accuracy of 
the present result is twice better than that of Ref. 113 811 and about two orders of magnitude higher 
than the experimental one. 



IV. CONCLUSION 



In this paper we have presented ab initio QED evaluations of the 2pi/2-2s and 2j)3/2-2s tran- 
sition energies in lithiumlike scandium, where the most accurate experimental data for middle-Z 
lithiumlike ions have been achieved. All presently available contributions to the transition energies 
are collected. Except for the one-electron two-loop correction, all other terms up to the two-photon 
level are treated within the framework of bound-state QED to all orders in aZ. The third- and 
higher-order interelectronic-interaction effects are accounted for within the Breit approximation 
using large-scale CI-DFS calculations. The relativistic recoil corrections are evaluated as well. 
As a result, the total theoretical accuracy is improved by a factor of 5 compared to the previous 
calculations. 

The ground-state hyperfine splitting of lithiumlike scandium has been calculated. The 
interelectronic-interaction correction to the first order in 1/Z is evaluated within the framework 
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of QED. The higher-order electron-correlation effects are calculated using the large-scale CI-DFS 
method. The one-loop radiative corrections are calculated with an effective screening potential. 
The theoretical value of the hyperfine splitting is improved in comparison with the previous results. 
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